The sign problem across the QCD phase transition 



Jens O. Andersen"^'', Lars T. Kyllingstad'', and Kim Splittorff^ 



c^ : 
o : 
o ■ 

(N ■ 
D ■ 

^ ; 
in : 



in ■ 
I ' 

Oh' 
^ . 



G\ 
O 
0^ 
O 



"'Niels Bohr International Academy, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 
Department of Physics, Norwegian University of Science and Technology, H0gskoleringen 5, N-7491 Trondheim, Norway 
Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 
(Dated: September 15, 2009) 

The average phase factor of the QCD fermion determinant signals the strength of the QCD sign 
problem. We compute the average phase factor as a function of temperature and baryon chemical 
potential using a two-flavor NJL model. This allows us to study the strength of the sign problem 
at and above the chiral transition. It is discussed how the Ua{1) anomaly affects the sign problem. 
Finally, we study the interplay between the sign problem and the endpoint of the chiral transition. 
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I. INTRODUCTION 

Suppose that there was no sign problem in QCD at 
nonzero baryon chemical potential. Then, by means of 
lattice QCD, it would be possible directly to study the 
phase diagram of QCD. In particular, it would be of 
great interest to get a direct computation of the pre- 
ferred ground state of strongly interacting matter for 
low temperature and chemical potential of order Aqcd- 
This region is relevant for dense stars, and a number of 
theoretically motivated ground states has been proposed 
P, S, H, 3 — see Ref. Q for a review. Clearly, a reliable 
lattice QCD test of these theoretical predictions could 
lead to substantial progress. It would also be of rele- 
vance to study chemical potentials and temperatures in 
the range of 100-400 MeV where ultrarelativistic heavy- 
ion collisions are conducted. 

In order to show that it is the sign problem which 
is the main obstacle for such first principle numerical 
investigations of the QCD phase diagram, let us for a 
moment consider phase-quenched QCD. In the phase- 
quenched version, one replaces the fermion determinant 
by its absolute value. Standard lattice techniques can 
therefore be applied. Studies of this theory on the lattice 
show a rich phase diagram [1, 0, [l] • For two flavors the 
phase-quenched theory is realized by coupling the chemi- 
cal potential to the third cornponent of the isospin rather 
than to the quark number [9|, and the phases obtained 
numerically can be understood physically. For example, 
the phase quenched lattice data show that a second-order 
phase transition takes place at zero temperature once the 
isospin chemical potential reaches half of the pion mass 
[1, 0, 111 ■ This is understood to be the transition into a 
phase with a Bose condensate of pions. Besides offering 
direct insights into an alternative direction of the QCD 
phase diagram, the phase quenched studies also demon- 
strate that without the sign problem it would be possible 
also to study nonzero quark chemical potential by means 
of lattice QCD. 

While a direct solution of the QCD sign problem eludes 
us at present, progress has been made to circumvent the 
sign problem 0, lil, E, M M M M [13, M M M 



|2l|, 122, l23|, l2J, l25|, l26|. The success of these approaches 
is linked to the degree of cancellations which take place 
due to the sign problem. One way to measure the degree 
of such cancellations is by means of the average of the 
phase factor 
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If the average of the phase factor is tiny then severe 
cancellations take place in the path integral defining the 
QCD partition function and the validity of lattice com- 
putations should be examined carefully. Studies of the 
average phase factor are therefore important in order to 
understand lattice gauge computations which circumvent 
the sign problem. 

Chiral perturbation theory has been used to predict 
the behavior of the average phase factor [23, [H, l29| . 
These predictions have been confirmed by lattice QCD 
measurements. Of course, as the temperature approaches 
the critical one, chiral perturbation theory is no longer 
valid since the more massive modes arc no longer Boltz- 
mann suppressed. For example, it has been demonstrated 
in Ref. [3l| that the effect of the nucleons on the average 
phase factor can be tracked in lattice QCD simulations 
just below the critical temperature. 

The purpose of the present study is to examine the 
average phase factor also at and above the chiral phase 
transition, where we can not trust chiral perturbation 
theory. In order to carry out such a study we necessarily 
need to rely on a specific model. Here we have chosen the 
two-fiavor NJL model. This is perhaps the most widely 
used framework to model asp ects of the QCD phase di- 
agram [3^ . Issl . Isl . Issl . [H, [sj [111 . For a review, see Ref. 

Since the phase factor is what separates phase quen- 
ched QCD from full QCD, the expectation value of the 
phase factor in the phase quenched theory is simply 
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where the subscript l-f- 1 refers to the ordinary two-flavor 
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theory and 1 + 1* refers to the two-flavor phase quenched 
theory (since it has a flavor and a conjugate flavor). 

In the present paper we will evaluate this ratio within 
the NJL model at mean field level. The general structure 
of the saddle point approximation gives 
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where D, is the mean field free energy density and H is the 
Hessian matrix describing the fiuctuations at the saddle 
point. 

For small values of the chemical potential, the mean 
free energies in the 1 + 1 theory and the 1 + 1* theory 
are identical. In this case, the prefactor — the ratio of the 
Hessian determinants — determines the average phase fac- 
tor. For /i > 77i7r/2, the phase quenched theory enters the 
Bose condensed phase mentioned above. The two mean- 
field free energies are therefore different in this phase. 

This effect was first understood within chiral pertur- 
bation theory. It is also present when the average phase 
factor is evaluated within a chiral random matrix model 
[40| and also results from an exact computation in the 
e-regime of chiral perturbation theory [U, \^ . 

The Ua{^) axial anomaly affects the order of the chiral 
phase transition at zero chemical potential [i^ . A recent 
study [13] shows that the location of the endpoint of the 
first order chiral transition can be highly sensitive to the 
strength of the Ua{^) breaking. Here we compute the 
effects of the C/a(1) breaking on the average phase factor 
within the two-fiavor NJL model. 

Since the strength of the sign problem sets a practical 
upper limit for lattice QCD it is crucial to understand 
the value of the average phase factor at the endpoint. 
Here we discuss the interplay between the endpoint and 
the average phase factor. Perhaps somewhat surprisingly 
we show that the average phase factor grows (indicating 
a milder sign problem) as one approaches the endpoint. 
The reason is that the massless mode associated with the 
endpoint induces a pole in the prefactor of Eq. ([3]). 

The paper is organized as follows. In Sec. II, we briefiy 
discuss the two-fiavor NJL model and its symmetries. In 
Sec. HI, we study the average phase factor without axial 
symmetry breaking, and in Sec. IV, we investigate the 
effects of axial symmetry breaking. In Sec. V, we look at 
the interplay between the sign problem and the endpoint 
of the chiral transition, and in Sec. VI we summarize our 
results. 



II. THE NJL MODEL AND THE AVERAGE 
PHASE FACTOR 

In this section, we briefiy discuss the two-flavor NJL 
model and our conventions. Note in particular that we 
are including the prefactor in our mean-fleld calculation 
of the partition function, cf. Eq. (|28p. 

The fermionic part of the two-flavor QCD Lagrangian 
at flnite baryon chemical potential and flnitc isospin 



chemical potential is 

iJ>{D - m + + M/^37o)i/', 
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where B = diag(l/3, 1/3), I3 = T3/2, and D is the 
QCD fermion operator at zero chemical potentials. Since 
the Lagrangian is quadratic in the fermionic fields, the 
fcrmions can be integrated out exactly, yielding the 
fermion determinant 

det(L' + ?7i - /iB-B7o - Ai//37o) 
= det(D + jTi — ^-fo — (5^70 ) 

X det(Z) + m - /.(7o + (S/i7o) , (5) 

where the chemical potentials /i and 6fi are 



LL = ^-^ and 5ii — —. 
'^3 ^2 



(6) 



/i is thus the quark chemical potential. 

As mentioned in the introduction, phase-quenched 
QCD with two fiavors is realized by setting = and 
interpreting 5^ as the quark chemic al p otential instead. 

The NJL Lagrangian is given by [391 



where the various terms are 
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In the two-fiavor case, ip is an isospin doublet, i.e. 



(11) 



and m ~ dia,g{mu,md) is the quark mass matrix. In the 
following, we take niu = md = mo. Also, r = (ti, r2, ra), 
where the r^s are the Pauli matrices acting in isospin 
space. 

The Lagrangian ([7|) is invariant under global SU{Nc) 
and U{1)b symmetries. The latter corresponds to conser- 
vation of baryon number. The interaction terms Ci and 
£2 both have a a global SU{Nf)i, x SU{Nf)Yi, which is 
broken explicitly to a 5[/(A^/)l+r by a nonzero value 
of TTiQ. Finally, £1 has a U{1)a symmetry which is not 
shared by £2 . This means that we can use the parameter 
a to control the strength of the axial symmetry violation 
in the model, and when a = 0, the Lagrangian ^ is 
U{1)a symmetric. 

We introduce chemical potentials for baryon number 
density and isospin density by, respectively, adding the 
terms 



£b = y^Bi^l^Bip, 



(12) 
(13) 
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Using the definitions in Eq. ([6]), we can now write the 
Lagrangian as 



+G [iH'f + (^«7^T?A)2] 

+ (1 - 2a)G [i^r^ f + (^nV)'] 



(14) 



When a = 1/2, the terms in the scalar- isovector and 
pseudoscalar-isovector channels vanish, and one recovers 
the standard NJL Lagrangian. 

Next, we introduce four auxiliary scalar fields a and 
fli, and four pseudscalar fields rj and tt^, (i = 1,2,3), by 
adding the following terms to the Lagrangian: 



C 1- 

1 - 
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{f] - 2G^ij^^f + (a, - 2G^nijy 



(15) 



If we use the equation of motion for the auxiliary fields 
to eliminate them, we recover the original Lagrangian. 
The quartic interaction terms in the Lagrangian vanish 
and we obtain 



-ij^TiTTi - (1 - 2a){ij^r] + r^a;)] ip 



4G 



(cr^ + TTiTTi) 
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AG 



{if+a^a,). (16) 



Based on the last term in Eq. it is evident that we 
must require a < 1/2 for the theory to be stable, which 
means that the Uyi(l) symmetry is maximally violated 
for a 1/2. 

To allow for a chiral condensate (V'V') a-nd a charged 
pion condensate (ipij^Tiip), we introduce nonzero expec- 
tation values for the fields a and tti. Expanding around 
these values, we have that 



a-2G{^ij), 

TTi - 2G(V'i7ViV'). 



(17) 



For notational simplicity, we introduce the quantities 

M = mo-2G{^p^), (18) 
p = -2G{4)i-i^Ti^). 

We can now write the Lagrangian as 

C = [i7^9^ - A/- ^7^1/9 + ^^7" + '5Ai7°T-3 

-cr - ij^n-Ki - (1 - 2a){ij^T] + Titti)] ip 
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The pion condensate p breaks parity and isospin symme- 
try. 

Note that when a ^ 1/2 we should, in general, allow 
for nonequal (uu) and (dd) condensates, since (03) = 



~2G{{uu) — {dd)). However, in this paper we will only 
consider the special cases {p, ^ 0,6p = 0) and {p = 
0, 6p ^ 0), and in both these cases (uu) — (dd). 

After applying the shifts in Eq. (??), the Lagrangian 
is quadratic in the fermion fields, which means that they 
can be explicitly integrated out. Thus, the partition func- 
tion of the NJL model is 

» 3 3 
^NJL - / VaY[Vn,Vi^ l[Va, e"^""!--'^-'''''^!, (20) 

z=l j=l 

where the effective action is given by 

<Soff[o-, TTi,??, Oi] = -logdeti^T 
fP r 

+ TTiTTi + {I - 2a){if + a^ai) , 

(21) 

and the matrix K is defined as 

K = 7^P^ + niQ - p-/o - Sp-foT3 + a 

+i-izmii + (1 - 2a)(i7^7 + nai) . (22) 

The matrix K has components in Dirac, color, and 
isospin space, as well as spacetimc, and the determinant 
must be taken over all these spaces. In general this is 
highly nontrivial to do analytically, but it can be simpli- 
fied somewhat by defining a matrix K' by 



55 /^'-i * 



(23) 



Since det(7o) = 1, this matrix has the property that 
log(detiir') = log(detiir). Unlike however, K' has 
all pqS on the diagonal, so it can be written as 

K' =p^ + K'{p^ = Q). (24) 

Thus, the eigenvalues of K' (and K) can be written as 

K=Po+£i, (25) 

where are the eigenvalues of the simpler matrix 
K'{po = 0). The fermion part of the effective action 
is then 



log detK = ^ log(po + £»)■ 



(26) 



Since the original matrix K is the inverse fermion prop- 
agator, we immediately identify Si as the energies of the 
quasiparticles in the spectrum. 

The values of the condensates M and p are found by 
minimizing the effective potential in the (M, p) plane, i.e. 
by simultaneously solving the equations 
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(27) 



while ignoring the fluctuations in the fields. 
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Using the saddle-point approximation to evaluate the 
partition function at this minimum, we find that 



1 



VdctH 



-pvn 



(28) 



where the thermodynamic potential O is the effective ac- 
tion evaluated at the saddle point and divided by the 
four- volume PV . H is the Hessian matrix of ScS, defined 
as 
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(29) 



and also evaluated at the saddle point. This matrix is a 
measure of the curvature of the potential at the saddle 
point, i.e. the magnitude of the fluctuations of the fields 
around this point, and it is thus related to the masses of 
the effective excitations. 

Using the NJL model partition fmiction, we can now 
write Eq. ^ as 



(30) 



where the prefactor F is the square root of the ratio be- 
tween the Hessian determinants in the 1-1-1* and the 
1 H- 1 theories, i.e. 



_ det H{^i = 0,S^ ^ ^iq) 
y det H{fi ~ ^q^Sfi — 0)' 



(31) 



while Afl is the difference in thermodynamic potential 
between the two cases: 

Af2(/iq, T) = 51(/i = /iq, = 0) — il{fl = 0, 5/i = ^q). 

(32) 

Here, fiq is used to denote the quark chemical potential, 
so that fiq = /xb/3. 

In the numerical analysis, all momentum integrals have 
been regulated using a sharp ultraviolet cutoff of A = 
651 MeV. We have taken the bare quark mass and the 
coupling constant to be mp = 5.5 MeV and G = 2.12/A^, 
respectively, and set the number of colors to A^c = 3. 
With this parameter set, the model reproduces the pion 
mass TOtt = 139 MeV, the pion decay constant /^r — 94 
MeV, and a chiral condensate of {uu) = —(250 MeV)"^ in 
the vacuum. 

In the plots of the average phase factor (e^'^) we have 
set the three-volume V equal to L^, where L = 4/T. 
Thus, the full four-volume is 



pv 



64 

rp4 ■ 



(33) 



This choice is inspired by lattice QCD, and corresponds 
to a lattice size of 4 x 16^. 



III. THE IMPORTANCE OF THE PREFACTOR 

In Fig. [1] we show the chiral condensate as a func- 
tion of function of T/m.^ and ii/m.^. In Fig. [2] we show 
the chiral condensate as a function of T/niT^ and /Lt/m^. 
Comparing the two, we immediately see that the free en- 
ergies must differ for fi, Sfi > and low T. Indeed 
this is confirmed by computing directly the difference of 
the two free energies. This is shown in Fig. [31 This differ- 
ence is due to the formation of a pion condensate, p in the 
phase-quenched theory. In Fig. [4] we show a plot of the 
pion condensate as a function of T/m-^ and 5p,/m-^. As 
expected, the pion condensation sets in at 5p, ~ for 
T = 0. As the temperature is increased it takes a some- 
what larger chemical potential to trigger the formation 
of the pion condensate. As the pion condensate forms, 
the chiral condensate starts dropping to zero. 
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Figure 1: Tlie chiral condensate M in the 1-1-1 theory with 
a = as a function of T /rriTt and fi/m-,r- The curve marks the 
first order phase transition, while the dot marks the critical 
point. 

We conclude from this that the average phase factor is 
exponentially small when there is a pion condensate in 
the phase-quenched theory, cf. Eq. (j30p . However, when 
the pion condensate vanishes in the phase-quenched the- 
ory, the mean field NJL free energy is equal to that for 
the ordinary theory, i.e. Ari(/j,g,T) = if p = 0. The 
average phase factor, Eq. ((30)) . is then dominated by the 
prefactor F throughout an extended region of the (/ig, T) 
plane. 

The Hessian matrix measures the curvature of the po- 
tential at the minimum and can hence be interpreted as 
a product of quasiparticle masses. For low temperatures 
where the pions dominate the partition function, the pre- 
factor is simply [27j 



F = 1- 



2flq 



(34) 



This low-tcmpcraturc behavior of the prefactor is con- 
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Figure 2: The chiral condensate M in the 1 + 1* as a function 
of T /m-n and theory with a = 0. 
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Figure 3: The difference in the free energy densities (ASl), 
between the 1 + 1 theory and the 1 + 1* tlieory as a function 
of T/m^ and /i/m^ with a = 0. 



firmed by the NJL model computation: In Fig. Owe show 
the prefactor F as predicted in chiral perturbation the- 
ory (solid curve) and the prediction by the NJL model 
(crosses). 

The full plot of the average phase factor in the (jUg, T) 
plane is given in Fig. [6l We see that indeed the sign 
problem is exponentially bad in a region corresponding 
to the region of pion condensation in the phase quenched 
theory. Before reaching this the prefactor dominates and 
the average phase factor drops smoothly from one to zero. 

The full dominance of the prefactor outside the pion- 
condensed region is an artifact of the mean-field approx- 
imation. This has been demonstrated for low tempera- 
tures in chiral perturbation theory, where one-loop pion 
effects give rise to a nonzero value of Af2 also where the 



Figure 4: Tlie pion condensate p in the 1 + 1* theory as a 
function of T /rriTT and &p./m-n. with a = 0. 
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Figure 5: The prefactor F for T = 1 MeV and a = (crosses) , 
and the prediction 1 — (2/iq/m-n-)^ from Ref. [S^] (line). Note 
that the choice of T = 1 MeV instead of T = is purely for 
computational convenience. 

pion condensate vanishes [27j . By comparison to lattice 
QCD measurements of the phase factor, it is seen that 
both the prefactor and the one-loop contribution to Afi 
are relevant for the standard lattice volumes used [28t . 
Nevertheless, the overall trend for the average phase fac- 
tor also above the critical temperature is in agreement 
with lattice studies [isl ]. 



IV. THE EFFECT OF Ua{1) SYMMETRY 
BREAKING 

Instanton effects are responsible for the Ua{^) axial 
anomaly in QCD and is known [ist to have a strong in- 
fiuence on order of the phase transition at zero chemi- 
cal potential. In Ref. [43| it was demonstrated that the 
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Figure 6: The average phase factor (e'^'^) as a function of /i^m^ and T/m-^ for a = 0. 



strength of the ?7a(1) axial anomaly is also important 
for nonzero in the three-flavor NJL model. Here we 
consider the effects of the anomaly on the pha se factor 
within the two-flavor NJL model. As in Ref. [4J| we work 
with a jUg-dependent a, 



a 



1 2/2 



(35) 



so that the Ua{^) symmetry is restored at /x^ ^ /ig. The 
scale /xq introduces a new parameter in the NJL model. 
The form ([35|) is motivated by the Gaussian suppression 
of the instanton density due to Debye screening at large 
values of the chemical potential. 

The effects of a on the phase diagram of the two-flavor 
NJL model at nonzero ^ and 5^ have been discussed in 
R-cf. [s^. There, the focus was to demonstrate that the 
splitting of the {uu) and {dd) transitions induced by 5^ is 
wiped out by the breaking of Ua{^) symmetry. Here the 
setup is somewhat different since wc do not have /i and 
(5/i nonzero at the same time; recall that we are looking 
at the ratio Zi+i/Zi+i* . This implies that {uu) = {dd). 
Thus there is no effects of a in the free energies of the 
mean-field two- flavor NJL model considered here, and the 
dependence of the average phase factor on a is entirely 
due to the a-dependence of the prefactor F. 

In Fig. [71 we have plotted the average phase factor 
for four different temperatures, both with and without 
axial symmetry breaking. For the plot we have chosen 
/io = nriT, . We observe that there is no effect of a for low 
temperature whereas the Ua{^) breaking increases the 
average phase factor slightly for higher temperatures. 

From a Ginzburg-Landau perspective, it is easy to see 
that the effect of J7a(1) breaking on the phase diagram in 
the two-flavor theory is small, since the Ua{^) breaking 
only gives a quadratic term in the order parameter. This 
term is already present at a = 1/2. For Nf = 3 the 
Ua{X) breaking induces a cubic term in the Ginzburg- 
Landau free energy — see for example Ref. [i^ . This term 
potentially has a much larger effect on the position of the 




Figure 7: The average phase factor (e^*^) for T = m-^ (lower), 
T — 2m^, T = SrriTr, and T — 4m^ (upper), when a = (solid 

_ 2 / 2 

curves) and a — -^e ^q'"^-^ (dashed curves). 



endpoint. Indeed this is what was observed for the three- 
flavor NJL model in Ref. SI. 



V. THE ENDPOINT AND THE SIGN 
PROBLEM 

The endpoint of the first order chiral phase transition 
has been studied intensively in recent years. One con- 
clusion which has emerged is that the location of the 
endpoint in the {fiq,T) plane is extremely sensitive to 
the value of the quark masses and the fate of the ?7a(1) 
breaking. In order for lattice QCD to be able to reach 
the location of the endpoint in the {fj,q,T) plane, be- 
fore the numerical noise due to the sign problem wipes 
out the signal, the endpoint must be located outside 
the region where the average phase factor is exponen- 
tially suppressed. Whether or not this is the case is un- 
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Figure 8: The prefactor _F as a function of [iqlrxiT^ and T jm-^ for a = 0. Note the spike at the critical point, which extends far 
above the maximum value of F . 
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Figure 9: Contour plot of the prefactor F around the critical 
endpoint, with a = 0. The solid line is the first-order phase 
transition, and the dot marks the endpoint. The dotted lines 
are contours of equal F, each separated by a power of e. 



clear. Both the position of the endpoint as well as the 
region where the average phase factor is exponentially 
suppressed depend on the quark masses. In chiral Ran- 
dom Matrix Models [iO, and in the two-flavor NJL 
model, the location of the endpoint is always inside the 
exponentially suppressed region. The endpoints found by 
multi-parameter reweighting techniques in both Refs. [l9| 
and [20| are located where the exponential suppression 
sets in [i^. Here we discuss how the sign problem is 
correlated with the critical endpoint. 

If the endpoint is located in the region where the aver- 
age phase factor is exponentially suppressed it is not acci- 
ble by current lattice methods. If, however, the endpoint 
is located at small values of jiq, i.e. yUendpoint < ™ir/2, 
then the story is somewhat different. Two scenarios 
are possible depending on what happens in the phasc- 



qucnchcd theory: 

1 ) The phase-quenched theory has an endpoint at ex- 
actly the same location as the full theory. In this case the 
average phase factor will not show any particular sign of 
the endpoint. 

2) The endpoint in the phase-quenched theory is not 
at the same location as in the full theory. In this case the 
average phase factor can grow as one approaches the crit- 
ical endpoint. The reason is simple: just like the massless 
mode in the pion-condenscd theory drives the prefactor 
F to zero, see. Fig. [51 the massless a at the endpoint 
will induce a pole in the prefactor since the Hessian of 
the 1-1-1 theory is in the denominator, cf. Eq. ([3]). We 
can illustrate this by looking at the prefactor, F, in the 
two-flavor NJL model. In Fig. [51 we have shown the 
prefactor as a function of /^<j/m^ and T/m^r. We see 
that at the endpoint in the 1-1-1 theory the massless 
a has produced a peak in the exponential prefactor. In 
this case the peak is suppressed by the exponential factor 
even with a small volume. If, however, the endpoint is 
located outside the exponentially suppressed region, the 
peak due to the massless a may survive to much larger 
values of V. This is because the free energy difference 
Af2 can be much smaller here, since a cusp in fii+i does 
not necessarily induce a large Afl. Of course, it must 
always be true that < (e^*^) < 1, so the pole and the 
exponential factor must always be considered together to 
ensure this. 

As mentioned above it is unfortunately not possible to 
fully illustrate this most interesting scenario in the two- 
flavor NJL model used here. However, it is quite possible 
that it could be realized in the three-flavor NJL model 
with J7a(1) axial symmetry breaking. 

One of the characteristic features of the critical point 
is the divergence of the correlation length ^. Since the 
correlation length is inversely proportional to the sigma 
mass, while the Hessian matrix is roughly proportional to 
the squares of the masses, we expect that as wc approach 
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the endpoint, 

F^^. (36) 

In Fig. [9] we have shown a contour plot of the prefactor 
F in the region immediately surrounding the endpoint. 
A comparison of this figure to Fig. 3 of Ref . |471] , which 
is a similar plot of the correlation length ^, seems to con- 
firm Ec^. ((36|) . A notable feature of this plot is the fact 
that the prefactor — specifically, the contribution from the 
prefactor to the Zi+i partition function — is discontinu- 
ous across the phase transition. This seems to indicate 
that there is a zeroth-order phase transition in the theory. 
Fortunately, this discontinuity is entirely due to correc- 
tions to the partition function of order 1/V, so there is 
no inconsistency in the infinite- volume limit. 

VI. CONCLUSIONS 

We have presented the first study of the average phase 
factor of the complex QCD fermion determinant within 
a two-flavor NJL model. This study extends previous 
studies in several directions: We have considered tem- 
peratures not accessible to chiral perturbation theory, we 
have considered the effects of the Ua{^) anomaly on the 
sign problem and we have discussed the interplay be- 
tween the sign problem and the critical endpoint of the 
chiral phase transition. The NJL model certainly does 
not capture all effects of QCD and one should not take 
the precise numerical values given here as firm predic- 
tions. However, several general conclusions can be drawn 
from this study of the NJL model. 

First, we have seen that the exponential prefactor in 
the mean-field approximation dominates for a large range 
of /iq when the temperature is above the critical one. 
This is in complete agreement with the result obtained 



from a mean field random matrix approach [40|. For 
lattice QCD, this means that one should expect that the 
masses of the effective excitations have a large effect on 
the strength of the sign problem. 

Second, the average phase factor depends on the J7a(1) 
axial anomaly when the temperature is close to or above 
the critical one. The present study indicates that the 
sign problem may become milder when UAi^) symmetry 
breaking is taken into account. 

Third, we have discussed how the massless mode at 
the endpoint can drive the average phase factor to larger 
values. Since this indicates a milder sign problem at the 
endpoint, this observation is highly relevant for lattice 
investigations of the endpoint. Of course the free ener- 
gies in the phase-quenched and the full theory must be 
slightly different if the latter is to have an endpoint while 
the former does not. However, this difference need not 
be numerically large (as it is in the pion condensed re- 
gion) and the exponential suppression of the prefactor 
thus may only set in at large values of V. Unfortunately 
is it not possible to study this effect fully within the two- 
flavor NJL model. 

The order of the QCD phase transition depends on 
the number of flavors as well their masses [i^ . It would 
be therefore very interesting to extend this study to the 
three-flavor case where it also known that the effects of 
the Ua{^) axial anomaly are greater. 
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